Ground Return Impedance : Underground Wire 
with Earth Return 

By JOHN R. CARSON 

Synopsis: In certain transmission problems principally those relating to 
induction and interference phenomena, it is necessary to know the trans- 
mission characteristics of a circuit composed of an underground wire with 
earth return. These can be evaluated by well known engineering formulas 
provided the ground return impedance is known. The present paper gives 
the mathematical solution of this problem and shows that the ground return 
impedance is substantially independent of the depth of the wire below the 
surface. 

THE object of this note is to give the solution of a problem of 
considerable interest and practical importance which does not 
appear to have been solved heretofore; this is the "ground return" 
impedance, per unit length, of a circuit composed of an underground 
wire or cable with earth return. 

The physical system and the problem may be more explicitly 
described and explained as follows: An underground wire or cable 
parallel to and at depth // below the surface of the ground is surrounded 
by a concentric dielectric cylinder of external radius a. The earth 
then forms the return path for currents flowing in the wire. The 
ground return impedance Z a is then defined as the ratio of the mean 
axial electric intensity at the external surface of the dielectric sheath 
to the current flowing in the wire. 

When the earth extends indefinitely in all directions about the 

wire so that circular symmetry obtains, the problem is quite simple, 

and the formula for the ground return impedance, denoted in this 

case by Z °, is well known. In practice, however, we are interested 

principally in the case where the wire is close to the surface of the 

earth, so that the distribution of return current in the ground is 

anything but symmetrical. For this case the formula for the ground 

return impedance, which it is the object of this note to state and 

discuss, is 

Z = (1 + c)Z a \ (1) 

Here the correction term c, which takes care of the departure from 
circular symmetry, is given by 



In formula (2), 



c = — ^ ft 2 '* " • h= d »- ( 2 ) 
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a' = dija, 

h' = h^fa, 

a = Aw'Kco where X is the conductivity of the ground in 
elm. c.g.s. units, co is 2t times the frequency, 



i = a/- 1. 

K (a'i->[i) is the Bessel Function of the second kind; it is equal to 

^ Ho li) (a'i-\^) where iJo (l) is the Hankel function as defined by 

Jahnke u. Emde in their Funktionentafeln. Denoted by ker a' + 
i kei a' the function K (a'i-\R,) has been computed and tabulated by 
the British Association. The only restriction on formula (2) is that 
the radius a is supposed small compared with the depth h. 

Now the ground conductivity X lies between 10 -14 and 10~ 12 , while 
the depth h will not in practice exceed a few meters (h < 10 3 ). Under 
such circumstances, at ordinary frequencies, h' will be exceedingly 
small compared with unity, and a' still smaller. Consequently in 
evaluating the infinite integral in (2), it is permissible to take e-^'VJS+i 
as unity, since for \x > 2, the rest of the integrand converges as i/4n 3 . 

Now we have 



j o v^r+m vf+i-U, r^-^ + v^rij 



^ = 2 



and hence c of formula (2) becomes 



1 

C = iz 



2 K (a'i-yli) 

Furthermore since a' by hypothesis is very small compared with 
unity, we can replace K by its limiting form for vanishingly small 
arguments which is approximately 

log (I/O- 

We thus get, finally, the approximate formula, valid for most practical 
applications, 

z - = { i+ ibiW)} z '°- (3) 

The interesting and somewhat surprising feature of this formula is 
that the value of the correction term 1/2 log (l/a') likely to occur in 
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practical applications amounts at most to 0.05 to 0.10. On the other 
hand, with the wire close to the surface of the ground, the conducting 
area of the ground return path is just one half the area available when 
the ground extends indefinitely in all directions and the return im- 
pedance is Z g °. In other words, the departure from circular symmetry 
means only a very small increase in the ground return impedance. 
In fact this increase is so small and the ground conductivity actually 
so variable, the correction is hardly justified by the precision of the 
data, so that, in most engineering applications, we may take Z g as 
equal to Z g ° with an error probably less than that involved in other 
factors, and lack of precision in data. 

Derivation of Preceding Formulas 

The derivation of the preceding formulas is not without interest. 
Since, however, this derivation is, in general, an adaptation of the 
methods employed in my paper 'Wave Propagation in Overhead 
Wires with Ground Return' (B. S. T. J., Oct., 1926) it will be out- 
lined rather than given in detail. 

Take the axis of the wire as the origin and Y as the vertical axis; 
then the surface of the ground is the plane y = h. Let a unit current 
flow in the wire and take the axis of the wire as the Z axis. In the 
ground (p = Vx 2 + / S a) the axial electric intensity will be written 

E = Z {. .- K^p'Hi) + E' = E° + £', (4) 

where p' = VaVx 2 + y 2 and K is the Bessel function of the second 
kind, related to the Hankel function by the equation 

The first term on the right hand side of (4) represents the circularly 
symmetrical distribution which would alone exist if the surface of 
the ground were removed to an infinite distance, while E' is a secondary 
distribution due to reflection at the surface of the earth (y = h). 
Inspection of equation (4) shows that when p = a, E is the required 
return impedance Z. 
Strictly speaking E° should be written as 

E° = ^-=- {K (p'i-\(i) + hK^p'i^cos 6 

+ faK 2 (p'i^i) cos 20 + • • • } , 
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the harmonic terms representing secondary reflection at the surface 
of the dielectric cylinder (p = a). If a is made sufficiently small, 
however, the harmonic terms become negligible. In view of this 
fact, the large amount of tedious additional analysis required, if the 
harmonic terms are retained, is not believed to be justified by the 
practical applications contemplated. E° will therefore be taken as in 
formula (4). 

The secondary electric intensity E' can always be written as the 

Fourier integral 

/•» 

E' = I F(fi)e"' ,J '^+ i cos x' n dn, ^ y < h, (6) 
Jo 

where x' = xya, y' = ysa, and the Fourier function F is to be 
determined. For the formulation of the boundary conditions at 
y = h we also require the expansion of K (p'i\t) as a Fourier integral; 
the required expansion is * 

K ( P 'i4i) = f f —— . e-^^+'cosx'n dn, P > 0. (7) 
Jo Vm + i 

In the dielectric, the magnetic forces H x , H y are taken as 

<j>(fi)e~ vl1 cos xtx dp 

r y^h. (8) 

II V = — I 0( / u)e -1 "' sin xn dfi 
Jo 

In the ground, on the other hand, we have 

io>H x = -4-E, 
dy 

a (9) 

io)II v = —E. 
ox 

In order to satisfy the boundary conditions at the plane y = h, we 
equate H x as given by (8) and (9), and H y as given by (8) and (9). 
The explicit formulas for H x and H y are derived from (9) by substi- 
tuting the Fourier integral for Ko, as given by (7), in (4) and differ- 
entiating as indicated. 

The two equations resulting from equating the two expressions for 
H x and the two expressions for H v can be solved for the Fourier function 

* See note at end of this paper. 
7 
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F(n). With this determined the required impedance Z is simply 
by (4) 

z = z g <> + v - o (io) 

on the assumption that a' = a4a is quite small compared to unity. 
This gives formula (1). 

Note: The expansion (7), which is believed to be novel, was derived 
by a limiting process rather too long and unsuitable for inclusion in 
this paper. It and the following additional expansions are quite 
useful in certain problems on wave propagation. 

cos 6-K^p'iJi) = - -fi f er*J*& cos x'n dp., 

Jo 

in O-K^p'i-yJt) = - <l H -7^-= e-^^+* sin x' p. dp., 
Jo -V/r + * 



sin 
where 



cos = y/p = y/^Jx 2 + y 2 , 



sin 6 = x/p = xj-^x 2 + y 2 . 



